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Abstract. We use the exterior product of double forms to re- formulate celebrated 
classical results of linear algebra about matrices namely Cayley-Hamilton theorem, 
Laplace expansion of the determinant, Newton identities and Jacobi's formula for the 
determinant. This new formalism is then used to naturally generalize the previous 
results to double forms. 

In particular, we show that the Cayley-Hamilton theorem once applied to the second 
fundamental form of a hypersurface is equivalent to a linearized version of the Gauss- 
Bonnet theorem, and once its generalization is applied to the Riemann curvature tensor 
(seen as a (2, 2) double form) is an infinitisimal version of the general Gauss-Bonnet- 
Ghern theorem. In addition to that, the general Gayley-Hamilton theorems generate 
several universal curvature identities in the sense of Gilkey-Park-Sckigawa. The exten- 
sion of the classical Laplace expansion of the determinant to double forms is shown to 
lead to general Avez type formulas for all the Gauss-Bonnet curvatures. 
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1. Euclidean invariants for bilinear forms 

Let y4 be a square matrix with real entries of size n. Recall that the characterestic 

polynomial of A is given by 

(1) 

n 

Xa{X) = det(A-AJ) = (-l)"A"+(-l)"-isi(A)A"-i + ...+s„(A) = ^(-l)"-*s,(A)A"-\ 

i=0 

Where 51(^4) is the trace of A, Sn{A) is the determinant of A and the other coefficients 
Sk{A) are intermediate invariants of the matrix A that interpolate between the trace 
and the detreminant, we shall call them here the Sk invariants of A. 
Since similar matrices have the same characterestic polynomial, therefore they have as 
well the same Sk invariants. In particular, one can define these invariants in an invariant 
way for endomorphisms. 

Let (y, g) be an Euclidean vector space of finite dimension n and /i be a bilinear form 
on V . We denote by h the linear operator on V that coresponds to h via the inner 
product g. 
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We define the Sk invariants of tlie bilinear form h to be those of the hnear operator h. 
Note that in contrast with Sk(h) , the invariants Sk{h) depends on the inner product g. 
In order to make this dependence exphcit we shall use the exterior product of double 
forms [3]. 

Recall that for 1 < k < n, the exterior product = h...h, where the bilinear form h is 
seen here as a (1, 1) -double form, is a multilinear form determined by the determinant 
as follows: 

(2) h''{xi, ...,Xk,yi, ...,yk) = k\ det[h{xi,yj)]. 
In particular, we have 

h" 

(3) deth = * — . 

n\ 

Where * denotes the (double) Hodge star operator operating on double forms, see [3]. 
Consequently, using the binomial formula, the characterestic polynomial of h takes the 
form 

{h-XgT 



XhW =det(/i - \g) 



m 



^' i=0 

n 

= ^(-l)"~^s,(/i)A"-\ 

1=0 

We have therefore proved the following simple formula for all the invariants of h: 
Proposition 1.1. For each < k < n , the Sk invariant of h is given by 

(4) = 

Where * denotes the (double) Hodge star operator operating on double forms, and the 
products g"'~'' , h'' , g'^~^h'' are exterior products of double forms, where g and h are 
considered as (1, 1) -double forms. 

In particular, the trace and determinant of h are given by 

Si{h) = * < —h )■ and Sn{h) = *—. 

[ (n — 1)! J n! 

Let us note here that for any orthonormal basis (cj) of {V,g), Sk{h) coincides by defi- 
nition with the Sk invariant of the matrix {h{ei,ej)). In particular, is the deter- 
minant of the matrix {h{ei,ej)). More generally, we have the following lemma: 
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Lemma 1.2. Let (cj), i = l,2,...n, be an orthonormal basis of {V,g), h a bilinear 
form on V and k + r < n. Then for any subset {ii,i2, ■.■,ik+r} with k + r elements of 
{1, 2, n} we have 

(5) g^h"" (ei,,ei2, ...,ei,^^,ei,,ei2, ...,ei,^,,) = (A; + r)!s^ (/i(ei^, e^J) . 

Where Sr {h{ei^, eij) is the Sr invariant of the (k + r) x (k + r) matrix (/i(ei„, Cj^)) , for 
l<a,b<k + r,gis the inner product on V . The product g^N' is the exterior product 
of double forms as above. 

Proof. The characterestic equation of the {k + r) x {k + r) matrix {h{ei^, e^J) is given 
by 

det (/i(ei„, CiJ - XI) = det {{h - Xg){ei^,ei^)) . 
Consequently, formula ([2]) together with the binomial formula show that 

det (/i(ei^ , J - XI) = jj^^-^{h - Xg) [d, , e^^ , . . . , e^,^^ , e^, , e^^ , . . . , e^,^^ ) 

= \^(_l)^+''--i\k+r-i..k+r-ij^i / \ 

To complete the proof just recall that the invariant of the {k + r) x {k + r) 
matrix {h{ei^,ei^)) is the coefficient of (— 1)'^A^ in the characterestic polynomial 
det {h{ei^, e^J — XI) , that is, jj^^:p^_g^h^ (cj^, e^j, Cj^.^^, , e^j, 6,^.^^) . □ 

In what follows in this section we shall use this formalism of the Sk invariants to refor- 
mulate and then generalize celebrated classical identities for matrices namely Laplace 
expansion of the determinant, Cayley-Hamilton theorem, Jacobi's formula for the deter- 
minant and Newton identities. 

1.1. Adjugates of bilinear forms. Recall that the adjugate (called also adjoint or 
cofactor matrix) of a square matrix of size n is a new matrix formed by all the cofactors of 
the original matrix. Where the (zj) -cofactor of a matrix is (— l)*"'^-' times the determinant 
of the (n — 1 X n — 1) sub-matrix that is obtained by eliminating the z-th row and j-th 
column of the original matrix. 

We are going to do here the same transformation but on a bilinear form instead of a 
matrix. Precisely, let (cj) be an orthonormal basis of {y,g), for each pair of indexes 
we define, the (zj)-cofactor of the bilinear form h., denoted tn-i{h){ei, Cj) , to be 
(—I)*"''-' multiplied by the determinant of the (ra — 1 x n — 1) sub-matrix that is obtained 
after removing the z-th row and j-th column from the matrix {h{ei,ej)). Then we can 
use bilinearity to extend into a bilinear form defined on V . It is then natural to 
name the so obtained bilinear form tk{h) as the adjugate of h. 

The next proposition shows in particular that the bilinear form tn-i{h) is well defined 
(that is it does not depend on choice of the orthonormal basis) 
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Proposition 1.3. // * denotes the Hodge star operator operating on double forms then 

(6) ^"-^^^)= 

Proof. From the definition of tlie Hodge star operator once it is acting on double forms 
we have 

(7) * (/^"-^)(e.,e,) = ^-l^(/.-^)(*e.,*e,). 

The last expression is by formula ([2]) exactly equal to (—I)*"''-' multiplied by the deter- 
minant of the (ra — 1) X (n — 1) sub-matrix that is obtained from {h{ei, ej)) by removing 
the i-th row and j-th column. This completes the proof of the proposition. □ 

We define now higher adjugates of h as follows 

Definition 1.1. For 0<A;<n— Iwe define the k-th adjugate of h to be the bilinear 
form given by 

(8) tk(h) = ^ ^ * ig^-^-^h^^ . 

^ ' ^ ^ k\{n-l-k)\ ' 

Note that to(^) = fi' is the metric, tn-i{h) coincides with the above defined adjugate of 
h. The terminology "higher adjugate" is justified by the following fact 

The quantity — ^ — - — y {g^''^^^''h''^{*ei,*ej) can be seen, like in the case k = n — 1 

above, as (—I)*"'"-' multiplied by a higher determinant of the (n — 1) x (n — 1) sub-matrix 
that is obtained from {h{ei, Cj)) by removing the i-th row and j-th column. 

Remark. In view of formula [5] it is tempting to think that the "higher determinants" 
— 7 — - — ^(Q"'~^~'^h'')(*ei,*ei) coincide with the Sk invariant of the matrix. However it 

turns out that this not the case in general: The higher determinants are multilinear in 
the rows (or columns) of the matrix while the Sk invariants are not in general, as can 
one can see easily in the case of the trace si . 

1.2. Laplace expansions of tfie determinant and generalizations. The following 
proposition provides a Laplace type expansion for all the Sk invariants of an arbitrary 
bilinear form. 

Proposition 1.4. For each k, Q < k < n — 1 , we have 

(9) {k + l)sk+i{h) = {tk{h),h). 
Where (., .) denotes the standard inner product of bilinear forms. 
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Before we prove this proposition let us clarify its relation to the classical Laplace expan- 
sion (called also cofactor expansion) of the determinant. Remark that for k = n — 1 we 
have 

n 

nsn = {tn~.i{h),h) = ^ t„„i(/i)(ei, ej)h{ei, ej). 

Recall that by definition the factor t„_i(/i)(ej, Cj) is the usual (ij)-cofactor of the 
matrix {h{ei,ej)) and s„ is its determinant. We therefore recover the classical Laplace 
expansion of the determinant. 

Remark. Actually, Laplace expansion of the determinant is more refined. Precisely, for 
any i or j we have 

n n 

Sn = ^tn-iih){ei, ej)h{ei, Cj) = tn-i{h){ei, ej)h{ei, ej). 
i=i j=i 

Proof. Using basic properties of the exterior product of double forms and the generalized 
Hodge star operator, see [31 [S], it is straightforward that 

{t,{h),h) = *i{*t,{h)}h) = * = (k + l)s,^,{h). 

□ 



1.2. L Further Laplace expansions. Recall that the determinant of h is determined by 
/i". As the later expression can be written in several ways as a product h"-~'^h^ for each 
r, we therefore get different expansions for the determinant by blocks as follows: 

, , /i" /i'' (n-r)\r\, /i""'' h\ 

det h = *— = * — = ■ (*- -, — ) 

nl nl nl [n — r)l rl 

yr) ii<i2<...<ir ' ^ '' 

jl<j2<---<jr 

Where {ei, 62, e„} is an orthonormal basis of V, e(p) and e{a) are the signs of 
the permutations p = and a = {ji,.-.,jn) of (1, 2, tt.) . Recall that 

^ (eij, ...,ei,,ej,, ...,ejj (resp. (e,^^,, ei„, e^^^,, e^^) equals the determinant 

of the r X r sub-matrix (/i(ejj., Cj-J) for 1 < k,l < r (resp the determinant of the 
n — r X n — r sub-matrix (/i(ejj,, e^J) for p + 1 < k,l < n). Note that the second 
sub-matrix is just the co-matrix of the first sub-matrix, that is the sub-matrix obtained 
from the ambient matrix {h{ei,ej)) of size n after removing the rows and the 

columns ji, ...,jn. Let us mention here also that the original Laplace expansion is finer 
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than the previous expansion, precisely it says that 

det^= e(p)e(a)— (ei,,...,ei,,ej,,...,ejj— — — y (e,^^,, e,„, Cj^^,, e^J . 

*l<«2<---<«r 

Remark. One can write easily similar expansions for all the lower invariants of h. 
In fact, the product g^~^h^ can be written in different ways as g'^h'^ g'^~^~^h^~'^ for 
Q < q < k and Q < p < n — k . Precisely we have, for 1 < A; < n and for every < g < A; 
and < p < n — A; , the following expansion 

k\{n - k)\sk{h) = ^g'^-^h^ = * (^p/iV"''"^/^''"') = (^^/^^ = 

ii<i2<...<ip+q 

jl<j2<-<jp+q 

Where as above {ei, 62, 6^} is an orthonormal basis of V , e(p) and e(cr) are the signs 
of the permutations p = (zi, z„) and a = (ji, j„) of (1, 2, ri) . 

1.3. Girard-Newton identities. 

Proposition 1.5 (Girard-Newton identities). For Q < k < n — 1 , the trace of tk is 
given by 

(10) ctk{h) = {n-k)skih). 

Where c denotes the contraction map. 

Proof. Using basic properties of the exterior product of double forms and the the gen- 
eralized Hodge star operator, see [SIE], we immediately get 

ctk{h) = *g* tk{h) = * < — — \ = {n- k)sk{h). 



{n-k- 1)! k\ 



□ 



In order to explain why the previous formula coincides with the classical Girard-Newton 
identities, we shall use lemma fLSl which gives an explicit formula for the linear operator 
tr{h) that is associated to tr{h) via the metric g. Denote by h the linear operator 
associated to h and let pi =iiW . Lemma 11.81 shows that 

r 

il tr{h) = ^(-l)^S^„i(/i)pi. 

i=0 

Therefore we can reformulate the identity ( JTOj) as 

r 

^(-l)*Sr._i(/i)pi = (n - r)sr{h), 



i=0 
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or 

r 
i=l 

That are the celebrated classical Girard-Newton identities. 

Remark (Terminology). The transformations tj. are famous in the literature as New- 
ton's transformations. Up to the author's knowledge, it was Reilly [9] the first to call 
them as such (he treated only the case of diagonalizable matrices). He justifies this by 
the fact that they generate the classical Newton identities as above. With reference to 
the above discussion, the terminology adjugate (or cofactor) is in the author's opinion 
more appropriate. 

1.4. Higher adjugates and Laplace expansions. Let h he a. bilinear form on the n- 
dimensional Euclidean vector space {V, g) . We define , for < q < n and < r < n — g, 
the {r,q) adjugate of h denoted S(^r,q){,h) by 

(11) S^r,,){h) = — r- * {g-^-'h''). 

ql(n — q — ry. 

Note that S(i^g)(/i) = tg{h) is the adjugate of order g of /i as defined in subsection ILII 
and S(o,g)(/i) = Sq{h) is the Sg invariant of h. 

The higher adjugates S(^r;q){.h) of h satisfy similar properties like the usual adjugate of 
h which was discussed above. We list some of them in the following theorem 

Theorem 1.6. For any integers r and q such that < q < n and 1 < r < n — q we 

have 



General Laplace's expansion: 
iq + r)\ 



q+r{h) = {S(^r,q){h),h''). 



q\ 

• General Newton's identity: 

c {s(r,q){h)) ={n-q-r + l)s(r_i_g). 

Proof. To prove the general Laplace's expansion we again use the formulas of [3] to get 
a one line proof as follows 

\ql[n — q — ry. J \[q + ry.[n — q — ry. J ql ql 
In the same way we prove the general Newton's identity as follows 

(gU-q-rf^q \ / gU-q-r ^q 
~i7 ] =*9( -r, 
ql[n — q — ry. J \q\[n — q — r 

gn~q~r+lj^q 



q\{n — q — r + 1] 



{n-q-r + l) = {n-q-r + l)s(^r-i,q)- 
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□ 

Remark. (1) If the bilinear form h is diagonalizable, that is if there exists an or- 
thonormal basis (cj) of V such that h{ei,ej) = Xig{ei,ej) for all We call 
the real numbers Aj the eigenvalues of h. Without loss of generality we assume 
that Ai < A2 < ... < A„. 

It is not difficult to show that all the double forms S(^r,q){h) with r > 1 are then 
also diagonalizable in the sense that 

f 

Where ii < 12 < ... < ji < 32 < ■■■ < jr and the eigenvalues of are 
given by 

•^ni2...ir ~ ^ ] '^ii -^i2 • • • -^iq • 

ii<i2<...<iq 

{il J2,.-- Jq}n{il,j2v,«r} = </> 

(2) In some applications it is useful to find the determinant of the sum of two matrices 
or more generally the Sk invariant of the sum. Using double forms formalism as 
above one can prove easily in one line the following identity for arbitrary bilinear 
forms A and B once seen as (1,1) double forms and for 1 < A; < tt, 



i=0 



1.5. Jacobi's formula. 



Proposition 1.7 (A Jacobi's formula for the Sk invariants). Let h = h{t) be a one 
parameter family of bilinear forms on V then 

(12) Jt''^^^ = ^^'-'^^^'1I^- 

In particular, for k = n we recover the classical Jacobi's formula: 

(13) ±det{h) = {t^.,{h),^). 

Proof. Using basic properties of the exterior product of double forms and the generalized 
Hodge star operator, see O |5] , we get 

-Sk{h{t)) - I ■■ ' - - 



□ 



dt ^ ' dt\ {n-k)\k\J \{n-k)\k\dt^ 

' f g^'-^h^^^ \dh\ , dh^ 
\ \n-k)\{k-l)\ ) rftj=^'-^(^)'dl^- 
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Remark. If one allows the inner product 5^ on to vary as well, say g = g{t) , then 
at t = we have the following generalization of the previous formula: 

(15) " {tk-i{h), ^) + {tk - Skg, ^). 

The proof is similar to the above one. 

1.6. Cayley-Hamilton Theorem. It is now time to give a sense to the top tk{h) , that 
is tn{h) , where n is the dimension of the vector space V. To simplify the exposition 
we assume here in this section that h is a symmetric bilinear form. Recall that for 
l<A;<?i— 1 we have, see [S] 

It is then natural to define tn{h) to be 

(17) Uih) = Sn{h)g - . \„ c"-^/t". 

Recall that = n!s„(/i), and therefore h'^ = n\sn{h) *1 = Sn{h)g"' . Consequently we 
have 

(18) tn{h) = Sn{h)g - Sn{h) . \„ c"-V" = Sn{h)g - Sn{h)g = 0. 

[n — ly. 

It turns out that the previous simple algebraic relation t„(/i) = is equivalent to the 
celebrated Cayley-Hamilton Theorem as we will see below. 



Lemma 1.8. Let < r < n — 1, and tr{h) denotes the linear operator associated to 
trih) then 

r 

(19) tJ}^=Y,{-Vfsr-,{h)h\ 

Where h is the linear operator that is asociated to h . 



Proof. The previous lemma is a direct consequence of the following induction formula 
(20) tr{h) = Sr{h)I - htr-i{h). 

Where / is the identity operator. In order to prove the previous formula let us diag- 
onalize h and denote its eigenvalues by Ai, ■ ■ ■ , A„. Then the eigenvalues of tr{h) are 
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given by (they are indexed by i) 

AijAij.-.Aj^ =— ^ Aj^Aj2...Aj 



il<i2<---<ir 



(21) 



--— ^ AiiAi2...Ai, - — ^ AiAi2...Ai^ - . 



11^^12^... J^ir 



--Sr{h) - Xi 



i^{i2,...,ir} 



( 



i2<...<ir 
\i^{i2,...,ir} 



i2^...J^ir 

i<^{i2,---,ir} 



(h) -Xi Xi^.-.Xi 



Ai2...Ai 



□ 



It results from the above discussion that 

n 

(22) t^To =x;(-i)'^n-w^'- 

and therefore the vanishing of tn{h) is equivalent to the Cayley-Hamilton theorem: 

n 

(23) ^(-l)^s„„,(/i)/i' = 0. 

1=0 

1.6.1. Cayley-Hamilton theorem vs. Infinitisimal Gauss-Bonnet theorem. Let M be a 
compact smooth hypersurface of the Euclidean space of dimension 2n + 1 . Denote by B 
the second fundamental form of M and by Sk{B) its Sk invariant. For each k, < k < 
n, the first variation of the integral J^^ S2k{B)dvo\ is up to a multiplicative constant, 
the integral scalar product {t2k{B), B) , where t2k{B) is the adjugate transformation of 
B as above, see [S]. The later result can be seen as an integral Jacobi's formula. By 
Cayley-Hamilton theorem t2n{B) = and therefore the integral Jj^j S2n(-S)dvol does not 
depend on the geometry of the hypersurface. In fact, the previous integral is up to a 
multiplicative constant the Euler-Poincare characterestic by the Gauss-Bonnet theorem. 
In this sense, Cayley-Hamilton theorem is indeed an infinitisimal Gauss-Bonnet theorem. 



1.6.2. Further algebraic identities of Cayley-Hamilton type. We continue here to assume, 
as declared in the beginning of this subsection II. 6[ that /i is a symmetric bilinear form. 
In this case the higher adjugates S(^r,q){h) are then as well symmetric (r, r) double forms. 
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Formula (15) of [3] provides the following expansion for < r < n and 1 < q < n — r: 

(24) su.,)ih)= y , ^ a'+'-idh'^. 
^ ' ^'^'^ ' ^ i\q\(%\r-q)\ 

i=max{0,g— r} 

This new form of S(^r,q){h) allows to extend its definition to the higher values of q, namely 
for q equal to n — r + 1, ...,n. For instance in the top case q = n > 2 we define 

(25) .M)(/^) = Yl TIFX ^9'^''-''c'h-. 

i=n—r ^ ' 

Recall that /i" = Sn{h)g'^ and therefore c*/i" = Sn{h)j^^g'^~\ consequently we have 
for r > 1 

We have therefore proved that 

(26) S(r,n)(/i) = for all 1 < r < n. 

For r = 1 we recover tn{h) = that is the usual Cayley-Hamilton theorem for h. The 
next case is when q = n — 1 > 2 and 1 < r < n — 1 , here we set 

i=n-r-l ^ J \ ' 'J 

Next we are going to show that 

(28) S(r,n-i){h) = for all 2 < r < n - 1. 

In order to prove the above identities, first remark that h^~^ is a (n — 1, n — 1) double 
form on an n dimensionnal vector space, then using proposition 2.1 of [5] we can write 
f^n-i _ gn-2j^ some (1, 1) double form k on V . Consequently, using some identities 
from [3] we get for i < n — 2 the following 

d{h--') =d{(f^-'k) = *f/ * -^^{n - 2)1 

= *g' {-k + gck) (n - 2)! = (- * g'k + *g'+^ck) {n - 2)! 

g'^-'-^k g'^^'-^ck \ {i + l)\)n - 2)\ 



- t\{n - 2)! . + ,^ . + ^—^^—^g-'-'ck 

\ [n — t — 2)1 [n — I — 1)1 J (n — 2 — 1)! 

{n-2)W. „_,_2/, t{ck) 



{n — i — 2)\ \ n — i — V 

For i = n — 1 , we get 

^n-l^n _ 1)1)2 
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Consequently the formula above defining S(r,„_i)(/i) takes the form 

5 (/,) = ^^^+ y (-ly+n-l ^^^^) , 

i=n—l—r ^ ' ^ ' ^ ' \ 



a r — n ^ ^\\(n — i — '?\\ I n — ^ \ ^ 



^ (z + r — + l)!(r?, — z — 2)! / n — 1 \ ^ (z + r — n + l)!(n — z — 1)! / n 
Changing the index of both sums to j = z — n + 1 + r we immediately obtain 



It is then easy to check that the previous two sums are both zero for r > 2. 

In the same way we define S(r,n-j)(^) using formula (I24p and one can prove similarly, as 

in the cases where z = and i = \ above, the following general result 

Theorem 1.9 (A general Cayley-Hamilton theorem). For l<z + l<r<n — i we 
have 

S{r,n-i){h) = 0. 

In order to re-formulate the previous theorem in terms of operators, like in the usual 
form of Cayley-Hamilton theorem, one can use the following induction formula 

(29) 

•S(r-,(j) {h) (Cjj , . .. , , Cjj , . .. , Cj^) l.g) (/z) (Cjj j •• Cj^ ; ^^2 ) •• ^ir ) 

/l(6jj , 6jj ) 1) (^) (^ii ; • • • ) ) ^ii i • • • ; ^ir ) ■ 

Where {ci, e„} denotes an orthonormal basis diagonalizing h and zi 7^ Z2 7^ ••• 7^ V- 
Consequently, we obtain the following 

g 

(30) S(^,g)(/i) (eii,...,eij = ^(-l)^S(^,_i,g_fc)(/i)((ei2, e^J /i'=(eij, 

fc=0 

where an over bar is made to indicate the correponding associated linear operators. 
In particular, we have for all 1 < i ^ j < n 

<? 

(31) J^~^{e,,e,) = ^(-l)'=V^(e,)/i'=(e,). 

A:=0 

Taking q = n — 1, the previous general Cayley-Hamilton theorem asserts that 

Corollary 1.10. Let V be a real vector space of finite dimension n > 3 and let A 
be a (diagonalizable) linear operator on V . Denote by tr{A) = Yll=o(~^y^r~i{A)A'^ 
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the r-th adjugate of A (that is the k-th Newton transformation), then we have for all 
^ ^ ^ j ^ n the following Cayley- Hamilton theorem 

n-l 

(32) 5^(-l)%_,(A)(e,)A^(e,) = 0. 

fc=0 

Where {ci, 62, e„} denotes an orthonormal basis diagonalizing A. 



1.7. Skew-symmetric bilinear forms and the PfafHan. Let us suppose in this sub- 
section that the bilinear form h is skew-symmetric. Then h can be seen either as a 
2-form or as a (1, 1) double form like what we did up to here in this paper. 
Suppose dimy = n = 2A; is even. We know already that *^ is the determinant of h 
once h is seen as a (1,1) double form. However, if we perform the same operations on 
the 2-form h and of course with the ordinary exterior product of forms and the usual 
Hodge star operator we obtain a new invariant Pf{h) called the Pfaffian of h. Precisely, 
we have 

(33) Fm = 

It turns out that Pf(/i) is a square root of the determinant of h that is (Pf(/i))^ = det h. 
This identity can be quickly justified as follows: 

We proceed by duality, in the case where /i is a 2-form, then ^ ® ^ is an {n, n) double 
form. Since the space of (n, n) double forms on V is 1-dimensionnal vector space then 
it is proportionnal to the double *^ where in the last expression h is seen as a (1,1) 
double form. It turns out that the previous two (n, n) double forms are equal. To show 
the desired equation it suffices to compare the image of the two double forms under the 
generalized Hodge star operator. From one hand we have 

On the other hand we have 




= det h. 



The Pfaffian satisfies several similar properties to those of the determinant and one may 
use the exetrior product of exterior forms (as we did here for double forms) to prove 
such properties. 

2. Euclidean invariants for Symmetric (2, 2) Double Forms 



In this section we are going to generalize the previous results to symmetric (2, 2) double 
forms that satisfy the first Bianchi identity. Recall that a (2,2) symmetric double form 
is a multilinear form with four arguments that is skew symmetric with respect to the 
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interchange of the first two arguments or the last two, and it is symmetric if we interchage 
the first two arguments with the last two. The Riemann curvature tensor is a typical 
example. 

During this section R denotes a symmetric (2, 2) double form on the n-dimensional 
Euclidean vector space (V, g) that satisfies the first Bianchi identity. 



2.1. The /i2fc invariants vs s^. invariants. For each k, < 2k < dimV = n, we 
define (by analogy to the invariants of the previous section) the h2k invariant of R 
to be 

(34) 

In particular ho = 1 and hn = *R^ in case n = 2k. 

In the case where R is the Riemann curvature tensor of a Riemannian manifold, h2k{R) 
is know as the 2A;-th Gauss-Bonnet curvature. 

Remark. Suppose n = 2A; is even and define the /i^-characterestic polynomial of R to 
be hn (^R — j • A direct computation shows that 

1=0 ^ ^ i=0 ^ ^ 

i=0 ^ ^ 

The /i„-characterestic polynomial of R is therefore a polynomial of degree k in X and 
its coefficients are all the h2i{R) invariants of i?, in particular these invariants coincide 
with the elementary symmetric functions in the roots of the /z„-characterestic equation 
of i?. It would be therefore interesting to identify these roots and to relate them in 
particular to the eigenvalues of the linear operator associated to R. 



2.2. Adjugates of (2, 2) double forms. Using the same procedure of cofactors as in 
the previous section we obtain several adjugates of R. Let us here examine the following 

(35) N2k{R) = „ and T^^iR) = *- ^ 



(n-2A;-2)! ' {n - 2k - 1)1' 

Note that N2k{R) is defined for 2 < 2A; < n — 2 and it is a (2,2) symmetric double 
form on V like R that satisfies the first Bianchi identity. On the other hand T2k{R) is 
defined for 2 < 2k < n — 1 and it is a symmetric bilinear form on V . 
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Theorem 4.1 of [3] provides explicit useful formulas for all the N2k{R) and T2k{R) as 
follows: 

''^^^ ~ (2A:-2)! ~ (2A:-1)!^ ^ 2(2A0!^ ' 

^''^^^ = ~(2k}\^~ i2k-iy: 

In particular, T2{R) = — cR is the celebrated Einstein Tensor. The higher T2k{R) 
are called Einstein-Lovelock tensors, see [6]. 

2.3. Laplace expansion of the h2k invariants and Avez formula. 

Theorem 2.1 (Laplace expansion of the /i2fc invariants). For A < 2k + 2 < n we have 

(37) h2k+2{R) = {N2kiR),R). 

Proof. Using the results of [3] one immediately has 

gn-2k-2j^k+l ^ ^2k+2j^k+l 



{N2kiR),R) 



{n- 2k -2)1} {2k + 2)\' 
This completes the proof. □ 

As a consequence we recover Avez's formula for the second Gauss-Bonnet curvature as 
follows: 

Corollary 2.2 (Avez's Formula). For n > 4, the second Gauss-Bonnet curvature equals 

h^{R) = |i?|2_|c/?|2 + l|c2/?|2. 

Proof. A direct application of formula |36] shows that 

N2{R) = R-{cR)g+''^g\ 



Consequently, 



h^{R) = {N2{R),R) = - {{cR)g,R) + {^g\R). 



To complete the proof just recall that the contraction map c is the adjoint of the mul- 
tiplication map by the g, see [3]. □ 



In the same way one can prove easily the following generalization of Avez's formula, see 
0, 
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Corollary 2.3. For 4 < 2k + 2 < n , the {2k + 2)-th Gauss-Bonnet curvature is deter- 
mined by the last three contractions of R'' as follows: 

^2^+2 = ((2;t-2)!' " \2k-l)V + 

2.4. Girard-Newton identities. 

Theorem 2.4. Let c denotes the contraction map then we have 

cN2k{R) = {n-2k- l)T2k and cT^kiR) = {n - 2k)h2k. 

Proof. Using the identity c* = *g from [3], one easily gets the desired formulas as 
follows: 

''^-f^) = ' * (n -2^-2). = * (n-2^-l)!(n-2t-2)! = <" " " ^'^^ 
Similarly, 

^"-^^-^i^^ g^-'^'^R'^ {n-2ky. 

'^^'^''^ = ' * (n-2fe-l)! = * (n-2fc)!(n-2fe-l)! = " 

□ 

2.5. Algebraic identities for (2, 2) double forms. 

2.5.1. The case of even dimensions. Suppose the dimension of the vector spave V is 
even n = 2k. We shall now give a sense the top T2k{R) ■ Using formula [361 we naturally 
set 

„n ok „n— 1 ok 

(38) UR) = '-^g-^-^ 

nl [n — 1)1 

Proposition 2.5. Let R be a symmetric (2,2) double form satisfying the first Bianchi 
identity on an Euclidean space of even dimension n then 

(39) Tn{R) = 0. 

Proof. Note first that if n = 2A; then /i„ = *R'' and therefore R^ = *hn. Where * is 
the Hodge star operator acting on double forms. For any < r < n we have 

gU-r 

[n — r)\ 

That is 

c'R'^ h 
r! (n — r)! 
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Next using the definition of T„ above we easily get tliat 

Tn = Kg - Kg = 0. 

□ 

Remark. In tliis case wliere n = 2k one could using formula fl36l) define Nn{R) as well 
by setting 

A direct adaptation of the previous proof shows that we have the algebraic identity 
iV„,(i?) = 0. 

2.5.2. The case of odd dimensions. Suppose now the dimension of the vector spave V 
is odd say n = 2k + 1. We shall now give a sense the top Nn-i{R) ■ Note that the top 
T„_i(i?) is well defined and need not vanish in general. Using formula 136) we naturally 
set for n > 3 

.r rr.^ d^-^R^ d'-'^R^ e-^R^ 2 

40 iVn-iOR = 7 ^ - 7 1^,9 + TTi Tv/- 

[n — 6)\ [n — 2)\ 2[n — \)\ 

Theorem 2.6. Let R he a symmetric (2,2) double form satisfying the first Bianchi 
identity on an Euclidean space of odd dimension ri > 3 then 

(41) iV„_i(i?)=0. 

Proof. Let n = 2/c + 1 > 3, note that i?^ is a (n — — 1) double form on an n 
dimensionnal vector space, then using proposition 2.1 of [8] we can write R'^ = g"^~^D 
for some (1,1) double form D on V^. Consequently, using some identities from |3] we 
get the following 



c 



= * c/"~^(n - 2)\{-D + gcD) = {n - 2)\ (- * g'^-'^D + *g'^~'^cD) 

={n-2)\{n-3)\(^gD+^-^^^^g' 

After contracting the previous identity twice we get 

^2fc-i^fc = (n - 2)!(n - 2)! {D + {n - 2){cD)g) , and c^'^R'' = {n - l)!(n - l)!cD. 

Consequently we have 

, , c'^-^R'' c'^-^R^ c'^-^R^ 2 

Nn-i[R) =-, TTT - 7 -r-g + -r-g 

[n — 6)\ [n — 2)\ 2(n— 1)! 

= {n - 2)\ (d + 'f'^^ g - D-{n- 2)cDg + {n - l)^g) g = 0. 



□ 
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Remark. In dimension n = 3 the previous theorem read 

N2{R) = R-{cR)g + ^g' = 0. 

In the context of Riemannian geometry where R represents the Riemann curvature 
tensor the previous identity is equivalent to the vanishing of the Weyl tensor in 3 di- 
mensions, in fact in this dimension N2{R) coincides with the Weyl tensor. 

2.5.3. Algebraic scalar identities for (2, 2) double forms. Suppose the dimension n of 
our vector space V is odd, say n = 2k + 1 and as above i? is a symmetric (2, 2) 
double form that satisfies the first Bianchi identity. . Corollary 12.31 allows one to define 
hn+i{R) = /i2fc+2(^), precisely we set 

We are going to show that h2k+2{R) as defined by the previous equation is zero. 

We proceed as in the proof of Theorem 12.61 and using the same notations of that proof 

we have 

„2fc-l pfc 

cR) = -{n-2y.{D + {n-2)cDg,cR), 



\2k-l 

r^k^k 2p 2p 

(_n_ =(n- l)\cD— 

^ {2k)\ ' 2' ^ > 2 ■ 

Taking the sum of the above three equation we immediately prove the vanishing of 
h2k+2{R) ■ Thus we have proved the following scalar identities 

Proposition 2.7. Let R be a symmetric (2, 2) double form satisfying the first Bianchi 
identity on an Euclidean vector space of odd dimension n = 2k + 1 > 3 then 

n-3j^fc c^'-'^Rk n^-^R'' c^R^ 

' R) - (t T^v cR) + (t TTT' ^) = 0- 



\n-?>)V ' \n-2)V ' \n-l 
In particular, for n = 3 we have 

{R, R) - {cR, cR) + ^ [c^Rf = 0. 

Remark. It the context of Riemannian geometry, where R is the Riemann curvature 
tensor (seen as a (2,2) double form), the previous scalar curvature identities coin- 
cide with Gilkey-Park-Sekigawa universal curvature identities [2] which are shown to 
be unique. Also the identities of Proposition 12.51 coincide with the symmetric 2-form 
valued universal curvature identities of [2] where they are also be shown to be unique. 
The higher algebraic identities, that are under study here in this paper, can be seen 
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then as symmetric double form valued universal curvature identities in the frame of 
Riemannian geometry. 

2.5.4. Higher algebraic identities for (2,2) double forms. Let k>l and < r < n — 2k 
and let R as above be a symmetric (2,2) double form on the ri- dimensional Euclidean 
vector space (V, g) that satisfies the first Bianchi identity. We define the (r, 2A;)-adjugate 
ot R, denoted /i(r,2fc)(-R) , by the following formula 

(42) ^"-)(«) = (;rd-7^ 

Note that /i(r-2fc)(-R) for r = (resp. r = 1 , r = 2) coincides with h2k{R) (resp. T2k{R) 
, N2k{R) ). 

Theorem 4.1 of [3] shows that 

(43) /.,...(«)= t 

i=max{0,2fc-r} ^ ' 

This last formula f H3|) allows us to define h(^r,2k){R) for higher r 's that is for r > n — 2k. 
The following theorem which provides general identities and generalize Proposition 12.51 
and Theorem 12.61 can be proved in the same way 

Theorem 2.8. Let R be a symmetric (2,2) double form satisfying the first Bianchi 
identity on an Euclidean vector space of dimension n . 

(1) If n = 2k is even then 

h{r,n-2i){R) = for 2i + l <r <n-2i. 

(2) If n = 2k + 1 is odd then 

h^r,n-2i-i)iR) =0 for 2i + 2 < r < n - 2i - 1. 

Remark that we recover Proposition 12.51 for n = 2k, r = l,i = and Theorem 12.61 is 
obtained for n = 2k+l,r = 2,i = 0. 



2.6. Jacobi's formula for double forms. 

Proposition 2.9 (Jacobi's formula). Let R = R{t) be a one parameter family of (2, 2) 
double forms then 

(44) ±h2t{R) = {kN2k-2{R),'^). 
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Proof. 

Jt^^^^> ~Jt V* (n-2fc)! ) ~ * 1, {n-2k)\ H 



□ 



Remark. If one allows the scalar product to vary as well say g = g{t) the previous 
formula takes the following form at t = 0: 

(46) ^Mk{R) = {kN,,_,iR), ^) + {T,, - h,,g, 

dt dt dt 

The proof is similar to the above one. 



2.7. Algebraic identities vs. infinitisimal Gauss-Bonnet theorem. Let {M,g) 
be a compact Riemannian manifold of dimension n = 2k. Denote by R its Riemann 
curvature tensor seen here as a (2,2) double form and, and let h2r{R) be the corre- 
sponding Gauss-Bonnet curvatures as above. For each r, < r < n, the gradient of the 
Riemannian functional Jj^^j /i2r(-R)dvol at g, once restricted to metrics of unit volume, is 
equal to T2r{R) , where T2r{R) is h2r{R) adjugate of R as above, it is known in geometry 
as the Einstein-Lovelock tensor, see [6]. The later result can be seen as an integral Ja- 
cobi's formula. Consequently, the algebraic identity Tn{R) = shows that the integral 
Jj^ hn{R)dvol does not depend on the metric g of M . In fact, the previous integral is 
up to a multiplicative constant the Euler-Poincare characterestic by the Gauss-Bonnet 
theorem. 

Here again, as in the situation of section ll.6.1l the linearized version of the Gauss-Bonnet 
theorem is an algebraic identity for (2, 2) double forms. 



2.8. PfafRan of 4-forms. Let a; be a 4-form on the Euclidean vector space V . Remark 

that u can be naturally considered as a symmetric (2, 2) double form. 

Suppose diml^ = n = Ak is a multiple of 4. By definition we have * u'^'^ = hn{ijj) 

is the hn invariant of u once it is considered as a (2, 2) double form. However, if we 
perform the same operations on the 4-form u and of course with the ordinary exterior 
product of forms and the usual Hodge star operator we obtain a new invariant Pf(a;), 
which we shall call the Pfaffian of uj . Precisely, we set 

(47) Pf(c.) = 

Using the same argument as in section [L71 it is plausible that (Pf(a;))^ = hn{u) . It 
would be interestiong to investigate the properties of this invariant for 4-forms. 
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3. Higher adjugate transformations and applications 



We shall now in this section generalize the previous results to higher symmetric {p,p)- 
double forms. 

Let a; be a symmetric -double form satisfying the first Bianchi identity, we define 
its adjugate transformation of order (r, pg) to be 



n-pq-r^q 

(48) hr,pq)i^) = TT- 

[n — pq — r)\ 

Where {] < r < n — pq. The result is a symmetric (r, r) -double form. We remark that 
for uj = h a (1,1) symmetric double form that is a symmetric bilinear form we recover 
the invariants of section 1 as follows 

h{o,i)(^) = Q-^qW, h(^i,q){h) = q\tq{h) and \r,q){h) = q\s(^r,q){h). 

Furthermore, for u = R a, symmetric (2, 2) double form, we recover the invariants of 
section 2 as follows 

/l(0,2,)(i?) = h2q{R), Wl,2q){R) = T^qiR) aud /l(2,29) (i?) = N2q{R). 

One can generalize without difficulties the results of the previous sections to this general 
setting. First using Theorem 4.2 of [3] one easily gets 



pq 



(49) hr,pq){uj) = J2 



.1 Y+pg 

_2 /-P9+ic*a;5. 



i\(r — pq + i)\ 

i=max{0,pg— r} 

As a first result we have the following Laplace type expansion: 

Theorem 3.1. Let u be a symmetric {p,p) -double form satisfying the first Bianchi 
identity on an n -dimensional Euclidean vector space V , Let q be a positive integer such 
that n > 2pq then 

(50) = {hpq,pq){u;),u;'') = ^-^{du^,du^). 

Proof. From one hand we have 

gn-2pq^q / gn~2pq^2q \ c^P^ (^uj'^g) 



(51) {h,,q,pq,{u),u^) = ^-J^;-^r^') - ■ V(n-2M)! ; - (2m)! ' 

On the other hand we have 

/ gn-2pq^q \ _ ^ (_l)r-+pg gT^r^q 

^ ' * \ {n-2pq)\ / H r!~' 

Consequently it is straightforward that 

„2pq(, ,2q\ PI ( 1 Nr+pg PI ( 1 Nr+pg 
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□ 

Taking u = R a. {2,2) double form we get 

Corollary 3.2 (General Avez Formula). Let R be a symmetric [2,2) -double form sat- 
isfying the first Bianchi identity on an n -dimensional Euclidean vector space V , Let q 
be a positive integer such that n > 4q then 

(54) /i4,(i^) = f;^(c'•i^^c'■i^''). 

In particular, 
(55) 

h, = \R\^ - \cR\^ + \\c^R\' and /ig = " |ci?T + 7|c'i?T - —\c'R'\' + ^k^i^T- 
4 4 36 576 

Remark. The previous corollary improves a similar result of [3] obtained for the case 
n = 4q. 

Corollary 3.3. Let h be a symmetric bilinear form on an n-dimensional Euclidean 
vector space V , Let q be a positive integer such that n > 2q then 

(56) s2,ih) = J2^-^n^{c''h^c^h'^). 
In particular, 

(57) S2{h) = -\h\^ + \ch\^ands4{h) = j/i^p - \ch^\^ + ^\c^h^\^. 

The following proposition can be seen as a generalization of the classical Newton iden- 
tities 

Proposition 3.4. For 1 < r < n — pq we have 

c (/i(,.p5)(a;)) = {n-pq-r + l)h^r-i,pq){^^) 



Proof. Using the identity c* = *g of [3], one easily gets the desired formulas as follows: 

C {hr,pq)i^)) = C * 7 77 = *7 77 = {u - pq - r + 1) /l(^„i,pg) (w) . 

\ y I — pq — rjl [n — pq — ry. 

□ 

We finish this section by establishing higher algebraic identities. First note that formula 
(149|) allows to define h(^r,2k) {^) for higher r 's that is in the cases where r > n — pk. The 
following theorem provides general algebraic identities and generalize Theorems 12 . 8 1 and 
12.61 and Proposition 12.51 It can be proved by imitating the proof of Theorem 12.61 
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Theorem 3.5. Let u be a symmetric {p,p) double form satisfying the first Bianchi 
identity on an Euclidean vector space of dimension n . Then 

h(r,pk-pi){^) =0, for n — pk + pi + 1 < r < pk — pi. 
In particular, we have 

(1) If n = pk is is a multiple of p then 

h(r,n-pi){^) =0 for pi + 1 < r < n — pi. 

(2) If n = pk + 1 is then 

h(r,n-pi-i){^) = for pi + 2 < r < n — pi — 1. 

Remark that we recover the results of Theorem 12.81 for j> = 2, Proposition 12.51 for 
n = 2k, r = 1, z = and Theorem 12.61 is obtained for n = 2k + l,r = 2,i = 0. 

4. Final remarks and open question 

4.1. Adjugate transformations vs. Gilkey's restriction map. Following [H [2] 
we define be the space of invariant local formulas for symmetric {p,p) double 

forms that satisfy the first Bianchi identity and that are homogeneous of degree n in 
the derivatives of the metric and which are defined in the category of m dimensional 
Riemannian manifolds. In particular = Im,n is the the space of scalar invariant 
local formulas and If^ „ is the space of symmetric 2-form valued invariants. 
Recall that the homogeneity of order n for u = u{g) ^ 1?^n i^ equivalent to 

for all c 7^ 0. 

The last property implies in particular that if oj{g) G then its full contraction 

cPu{g) G Xm,n and &-^u{g) G X^ „. 

The Gilkey's restriction map r : Xm,n — ^ Xm-i,n is closely related to the adjuagate 
transformations as we will explain below. 

Let uj{g) G X^^^, recall that the adjugate transformation of uj is h(i^pq){u) = 

*7 G If, „ . For a tangent vector we have 

g{n-l)~pg^q 

That is the restriction of the invariant formula Q = * (n~pq)] ~ ^{pqy. ^ ^m,n to the 
(n — 1) dimensional orthogonal complement of the vector v. 

More generally, h(^r,pq){^) is given by r successive applications of Gilkey's restriction 
map to the invariant formula G I,n,n • 
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4.2. Spectrum of adjugate transformations of (2,2) double forms. Let R be 

a symmetric (2, 2) double form defined over the n-dimensional Euclidean vector space 
{V,g). Denote by Ai, A2, Aat, where = ^^^~^'> ^ the eigenvalues of the linear operator 
A^V — )■ A'^V that is canonically associated to R. For each k with 1 < k < n — 2 , the 
exterior product g^R is a {k + 2,k + 2) symmetric double form and therefore has real 
eigenvalues once considered as an operator A^^^V — ?■ A^^^V . The eigenvalues of the 
operator are expected to be polynomials in the eigenvalues of R. The question here is 
to find explicit formulas for the eigenvalues of g'^R in terms of Ai, A2, Aat. 
In the case k = n — 2 the answer is trivial as we have in this case one single eigenvalue 
and it is equal to J2iLi • 

To motivate this question let us recall that the Weitzenbock transformation of order p, 
2<p<n — 2,of the double form R is given by |1] 

9cR g""' 

{p-2)V 

The positivity of the transformation N'p{R) has important consequences in Riemannian 
geometry via the celebrated Weitzenbock formula. 

The true question is then to determine the eigenvalues of N'p{R) in terms of the eigen- 
values of R. 

More generally, what are the eigenvalues of the exterior products g'^'R'^ and the contrac- 
tions d'R'^ in terms of Ai, A2, Aat? 

In particular, find explicit formulas for the invariants h2k{R) , for /c > 1, in terms of 
Ai, A2, Xn- 



4.3. Geometric applications of the general algebraic identities. We have seen 
in the sequel of the paper that the identity tn{h) = for a bilinear form leads to a 
linearized version of the Gauss-Bonnet theorem for compact hypersurfaces of the Eu- 
clidean space. Furthermore, for a (2,2) double form R, the identity Tn{R) = , for n 
even and in the context of Riemannian manifolds, is an infinitisimal form of the general 
Gauss-Bonnet-Chern theorem for compact manifolds. 

It is then natural to ask what differential or topological consequences can be drawn from 
the identity Nn-i{R) for a compact Riemannian manifold of odd dimension n? 
Note that for a 3-dimensionnal Riemannian manifold, the identity Nn-i{R) = for the 
Rieman curvature tensor R is equivalent to the vanishing of the Weyl tensor. 
The same question can be asked for the other different higher algebraic identities estab- 
lished here in this paper. 



4.4. A PfaflBan for 2k forms. Let a; be a 2A;-form on an Euclidean vector space (V, g) 
of finite dimension n. Remark that u can be naturally considered as a {k,k) double 
form. 

Suppose diml/ = n = 2kq is a multiple of 2k. By definition we have = ^{o,n)('^) 

is the /i(o,n) invariant of u once it is considered as a {k, k) double form. However, if we 
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perform the same operations on the 2 A; -form u and of course with the ordinary exterior 
product of forms and the usual Hodge star operator we obtain another invariant Pf (w) , 
which we shall call the Pfaffian of u . Precisely, we set 

(58) Pf(u) = *'^. 

Using the same simple argument as in section [L71 it is plausible that (Pf (w)) = /i(o,n.) {^) ■ 
It would be interestiong to investigate the properties of this new invariant for 2 A; -forms. 
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